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Abstract. We experimentally investigate the non-equilibrium steady-state distribution of the 
work done by an external force on a mesoscopic system with many coupled degrees of freedom: 
a colloidal crystal mechanically driven across a commensurate periodic light field. Since this 
system mimics the spatiotemporal dynamics of a crystalline surface moving on a corrugated 
substrate, our results show general properties of the work distribution for atomically flat surfaces 
undergoing friction. We address the role of several parameters which can influence the shape 
of the work distribution, e.g. the number of particles used to locally probe the properties 
of the system and the time interval to measure the work. We find that, when tuning the 
control parameters to induce particle depinning from the substrate, there is an abrupt change 
of the shape of the work distribution. While in the completely static and sliding friction 
regimes the work distribution is Gaussian, non-Gaussian tails show up due to the spatiotemporal 
heterogeneity of the particle dynamics during the transition between these two regimes. 
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1. Introduction 

A basic concept for a system driven in a non-equilibrium process by the change of some external 
parameters is that of work. For mesoscopic systems, e.g. colloidal particles or biomolecules, 
the work spent in such a process becomes a fluctuating quantity which follows a probability 
distribution of finite width mm- During the last two decades, various non-equilibrium work 
relations, above all the Jarzynski [3] and the Crooks relations [1], have been established and 
shown to restrict the shape of this probability distribution depending on the underlying specific 
features of both the system and the non-equilibrium process. From a more fundamental point of 
view, these relations refine the second law of thermodynamics at the mesoscopic scale. Further 
exact statements involving the applied work are rare, even though it has been demonstrated in 
the framework of stochastic thermodynamics mmm that the aforementioned relations as well 
as different ones for other thermodynamic quantities [SElEllgKin], namely entropy production 
and dissipated heat, can be derived from a broader perspective mm- 

Experimental tests of non-equilibrium work relations have been carried out for a variety 
of different systems; for single colloidal particles in time-dependent harmonic [121 fT3l fH] 
and non-harmonic potentials naiB], biomolecules in folding-unfolding assays dzi nHi [ig Eoj, 



Figure 1. (a) Schematic illustration of the sectional view of the experimental setup, (b) 

Snapshot of the equilibrium crystalline structure of the colloidal monolayer sX, B = 0.5 mT, light 
potential depth SOksT and v = 0. The blue bar represents 20/rm. See text for explanation. 


mechanical torsion pendulums coupled to a heat bath [21] . and charge transitions in electronic 
devices [221123] . Common to all experimental systems studied so far is that they only consist 
of a small number of degrees of freedom, whose internal interactions are irrelevant. From an 
experimental point of view the reason for this is straightforward: controlling all external forces 
acting on a system of many coupled degrees of freedom during a non-equilibrium process in a 
well-defined way, i. e., measuring the applied work, is a huge challenge [24] . 

On the other hand, a non-equilibrium phenomenon where the concept of work plays a 
prominent role is friction. For atomically flat surfaces sliding against each other, friction 
results from the interplay between externally applied forces and the nonlinear interaction of 
a large number of degrees of freedom making up extended contacts at the interface. Inspired by 
simplistic models [25] . a 2D system which has attracted much attention in recent years in the 
held of tribology in order to investigate in a controlled manner the spatiotemporal dynamics of 
crystalline surfaces under friction consists of a monolayer of interacting particles suspended in a 
huid and mechanically driven through a periodic potential. Since experimental realizations [26] 
and numerical simulations [271 [28] of this system have successfully shed light on the mechanisms 
behind friction, it represents also an appropriate model to investigate the statistical properties of 
the non-equilibrium work done by a well-controlled external force on a system composed of many 
interacting degrees of freedom. So far, this kind of analysis has only been numerically carried 
out to characterize plastic depinning of interacting particles within a stochastic thermodynamic 
context [29] . 

Here, we experimentally study the huctuations of the work done by an external force 
on a mesoscopic system with many coupled degrees of freedom: a crystalline monolayer of 
magnetically interacting colloidal particles moving on a periodic light held under commensurate 
conditions. The dynamics of this system mimics the transition from static to sliding friction, 
where a solid surface (the colloidal monolayer) is driven across a corrugated substrate (the 
periodic light potential) by an external force. We investigate the role of several parameters 
which can inhuence the shape of the work distribution, e.g. the number of particles used to 
locally probe the properties of the system and the integration time to measure the work. We 
hnd that, when tuning the control parameters to induce particle depinning from the substrate, 
there is an abrupt change of the shape of the work distribution. While in the completely static 
and sliding friction regime, the work distribution is Gaussian, non-Gaussian tails show up due to 
the spatiotemporal heterogeneity of the particle dynamics during the transition between these 
two regimes. Finally, we discuss the asymmetry of the work distribution within the context of 
the non-equilibrium fluctuation theorems. 
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Figure 2. (a) Trajectories of 50 neighboring particles measured over 200 s at constant potential 
depth Vmax = SO/cfiT and different applied forces /. From left to right: / = 12.4 fN, 18.6 fN and 
24.8 fN. (b) Trajectories of 50 particles measured over 200 s at constant applied force / = 12.4 fN 
and different potential depths Vmax- From left to right: Vmax = SOfe^T, 20kBT and 15 A;bT. 


2. Experimental description 

Parts of the experimental setup have been described in detail elsewhere [26] and will be discussed 
only briefly. A schematic illustration of the experimental setup is shown in Fig. (D^a). The 
system consists of a monolayer of superparamagnetic colloidal particles with 2r = 4.5/zm in 
diameter (Dyna M-450 Epoxy, Life Technologies), suspended in a 2.3g/1 sodium dodecyl sulfate 
aqueous solution and situated in a sample cell of 10 /zm height. The total number of particles 
forming the monolayer is ~ 5000. Using video microscopy, we simultaneously track the 
center of mass of approximately 500 particles in the full field of view at 3.3 frames per second 
with a spatial accuracy of 40 nm m- The viscous drag coefficient of the particles in the 
solvent under this confinement, measured from their equilibrium mean-square displacement, 
is 7 = 6.2 X 10“® kgs“^. A coupling between the particles is obtained by a static homogeneous 
magnetic field B = Be^ applied perpendicular to the sample plane. This field induces a repulsive 
dipole-dipole interaction U{d) = /zo(xi?)^/(47r(i^) with the magnetic constant, d the particle 
separation distance and x — 3.1 x 10~^^Am^ T“^ a constant, which allows to control the stiffness 
of the colloidal monolayer. By interference of three laser beams (A = 1064 nm) a light held with 
hexagonal symmetry inside the sample cell is generated, corresponding to a 2D periodic potential 
landscape V, whose maximum depth is Vmax = SOksT. The potential prohle is described by 
the function 
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Figure 3. Mean work done by the external force as a function on the number N of particle 
trajectories counted in the monolayer for two different frictional regimes (a) / = 24.8 fN and 
Vniax = ^OksT (stick-slip motion) and (b) / = 12.4 fN and Vmax = (complete sliding), 

computed over different time intervals: r = 0.6 s (o), 1.2 s (□), 2.4 s (o) , 4.8 s (<), 9.6 s (>), 
19.2 s (x), and 38.4 s (-I-). The dotted lines represent Eq. (I8|), with no fitting parameters. 


The lattice constant a and the depth Vmax can be tuned by the intensities and the angles of 
incidence of the laser beams m- Before a measurement is performed, the colloidal monolayer 
is allowed to equilibrate at room temperature T = 298 ± 0.5 K for at least one hour in presence 
of a magnetic held oi B = 0.5 mT. The resulting homogeneous crystalline state with hexagonal 
symmetry exhibits a lattice constant of approximately lO/iin with an interaction potential of 
C/(10/zm) ~ bksT, as shown in Fig. HKb). The lattice constant a of the light potential V is 
adjusted to the same value, i.e. we focus on commensurate conditions, in order to resolve the 
transition from static to sliding friction of the colloidal monolayer [26]. 

3. Dynamical response under applied force 

The whole sample cell is displaced in the x-direction with velocity v by use of a piezo table in 
order to move the particles across V, as sketched in Fig. [I](a). The movement allows to create 
a controlled homogeneous how v = ve^ and to drive the particles into non-equilibrium steady 
states (NESS). By tuning the value of v we can mimic the response of the colloidal crystal to the 
effective external force f = qv. The reason of this choice is that many tribological processes take 
place under these conditions [SUISSIIM]) where an external force is applied to a crystalline surface 
with the purpose of moving it against a substrate, and where the surrounding huid, in our case 
the solvent, only plays the passive role of a thermal bath to keep the system at hxed temperature. 
In order to induce a transition of the monolayer in response to f from static to sliding friction, 
we perform two different experimental protocols. In the hrst, we keep the depth of the light 
potential at its maximum value Vmax = 30k bT and then we displace the cell at different velocities 
to tune the magnitude of the applied force / = |f|. In Fig. [2Ka) we plot some trajectories of 50 
neighboring particles moving according to this protocol. At small velocities, corresponding to 
values of / much smaller than the maximum restoring force exerted by the light potential ([T|), 
max{—VF} = 87rVmax/(0a) « 34 fN, the particles remained pinned by the potential wells, as 





































shown in the left panel of Fig [2]^ a) for / = 12.4fN {v = 200nms“^). Note that, even when 
immobile in average, the position of each individual particle fluctuates due to the collision of 
the solvent molecules. As / approaches values comparable to or larger than 87rV'max/(9a)) he. 
when overcoming static friction, the particles are able to move across the potential barriers, thus 
resulting in collective motion. For example, at / = 18.6fN {v = 300nms“^), the monolayer is 
partially depinned and some of the particles start to move in the direction of f, as shown in 
the central panel of Fig. El^a). For / = 24.8fN [v = 400nms“^, right panel of Fig. [5]|^a)), all 
the particles in the monolayer are able to move in the direction of f. Under these conditions, 
the particles undergo stick-slip motion because their mobility is hindered by the presence of the 
potential barriers [261133| . resulting in a mean particle velocity 105nms“^ < v. In the second 
protocol, we fix the velocity of the sample at v = 200nms“^, i.e. at constant / = 12.4 fN, and 
we vary the depth Vmax of the light potential. The response of the particles to this protocol with 
decreasing values of Vmax is qualitatively similar to that observed when increasing / at constant 
Vmax and also has a transition from a static to a sliding friction regime. This effect can be 
observed in Fig. [2]^b) for potentials depths 30 ksT, 20 ksT and 15 ksT, at which the maximum 
restoring force 8TTVmax/{0a) has the values 34 fN, 23 fN and 17 fN, respectively. We point out 
that, although qualitatively similar, the particle dynamics resulting from these two protocols 
are not completely equivalent. Indeed, close inspection of the trajectories in Fig. [2]|^a) and EJ^b) 
shows that in the second protocol. Brownian motion is more significant and the particle mobility 
is higher because the local confinement created by the substrate is reduced when decreasing Vmax- 
Therefore, the second protocol reproduces the effect of changing the roughness of the substrate, 
which in turn results in higher particle velocities compared to the first protocol for the same 
values of the parameter 9/a/(87rl4iax)- For instance, at 9/a/(87rUmax) = 0.72, the mean particle 
velocity obtained by means of the first protocol is only 26% of the velocity v of the sample cell 
(right panel of Fig. [2](a)), whereas in the second protocol, it almost reaches free sliding at 93% 
of V (right panel of Fig. EKb)). 

4. Stochastic Thermodynamics of the monolayer under applied force 

We focus on the work done on the colloidal monolayer of N interacting particles driven across 
the corrugation potential V by an external force f, which is the common situation encountered 
in many tribological problems |32l [33l I34j . We fist present the equations of motion for our 
specific experimental protocol under flow v, which allows to mimic in a controlled manner the 
dynamics under applied force f . Then, we derive the corresponding stochastic-thermodynamic 
quantities of the latter tribological process. When the i-th particle (i = 1,...,^) moves at 
instantaneous position r* = {xi,yi) and velocity hj = {xi,yi) in presence of flow v, the viscous 
drag force relative to the flow is 7(rj — v). In our system, there is no actual external force but 
only conservative forces derived from the magnetically-induced repulsive interactions and the 
periodic light field. In addition, each particle is subject to the random thermal collisions of the 
solvent molecules. Therefore, the dynamics of the i-th particle is described by the Langevin 
equation 


liri-v) = -ViE + Ci, ( 2 ) 

where E is the total potential energy of the system, which includes the light potential V, the 
pair-interaction potential U of all the particles and the confining potential exerted by the sample 
cell, Vconf, which maintains the monolayer in a packed configuration and prevents the particles 
at the boundaries from escaping from the monolayer due to the repulsive interactions 

N . N 

E = J2 yir^) + 2 E E - r.l) + Vconf, 

i=l i=l jjti 


( 3 ) 


whereas the fast interactions with the surrounding solvent molecules are modeled by a Gaussian 
white noise of zero mean and correlations {t')) = 2'ykBT6{t — t')6ij. Because of the 

structure of Eq. ([2]), the dynamics of every particle in response to the flow and in absence of 
an external force [H] is equivalent to the dynamics in response to an external uniform force, 
f = 7 v, and without external flovJ^ 


yr* = f - ViE + f,i. (4) 

Therefore, hereafter we will focus only on Eq. Q in order to study the stochastic 
thermodynamics of a the monolayer under external constant force and without flow. 

In the context of stochastic thermodynamics, the first law for the potential energy variation 
along a single stochastic realization of the dynamics of the system can be written as HIM] 

N 

dE = ^ViE-dr,, 

= dW-dQ, (5) 

where dQ and dW are the heat dissipated into the solvent and the work done by f, respectively, 
and are given by 


N 

dQ = ^[f-V,E]-dri, 

2=1 

N 

dW = ^f-dr*. (6) 

2=1 

Then, from Eq. Q the work done by the uniform force f on the colloidal monolayer, normalized 
by ksT, over the time interval [0, r] reads 


Wr 


1 

/ 

ksT 


N 

^ f • i-j dt, 

2=1 
N 

^[a:i(r) -Xi(0)]. 
2=1 



(7) 


It should be noted that the expression of the work in Eq. ([7|) only involves the value of the 
force /, which can be tuned experimentally by means of u, and the instantaneous values of the 
x-coordinates of each particle, which are determined by videomicroscopy. Consequently, the 
work can be directly determined from the particles’ trajectories without the need to measure 
the pair interactions. Erom Eq. m, we can conclude that, regardless of the nature of the pair 
interactions, the mean value of the NESS work done by f over a time interval of duration r can 
be expressed as 


Hr = (Wt) 


Nf{x)T 

ksT 


( 8 ) 


where the brackets stand for an ensemble average over N particle trajectories and (x) = 
W is the drift particle velocity in response to f. We check that Eq. ([8]) is valid in 


^ This equivalence is only valid at sufficiently low Reynolds number, where the flow field around the particle is 
Stokesian and therefore the resulting drag force can be written as 'y{ii — v). In our experiments this assumption 
is fully justified because the Reynolds number is Re < 10““*. 





all the frictional regimes investigated in our experiments 0. For instance, it is trivially satisfied 
for static friction, where (x) = 0 yields fj-r = 0 because in average no mechanical work is done 
by f on the monolayer. On the other hand, for (x) > 0 the linearity of with respect to N 
and T predicted by Eq. ([8|) is also verified. For example, in Figs. Ha) and we plot for 
different integration times r the value of the mean work for stick-slip motion and free sliding, 
respectively, as a function of the number N of particle trajectories used in the computation of 
Wr- In this case, fir is determined by taking the average over all possible values of Wr at fixed 
N and r. We also plot as dotted lines the values of fir computed by means of Eq. dS]), where the 
NESS drift velocity (x) is independently determined from the particle dynamics. We observe 
that the agreement between both kinds of calculations is excellent. Therefore, from the validity 
of Eq. ([5]) we conclude the mean work mirrors the bulk frictional properties of the monolayer, 
namely a smooth transition from fir = 0 (static friction with zero mobility at small f) to fir <x f'^ 
(sliding friction with constant finite mobility at sufficiently large /) |26j . 


4.1. Non-interacting particles 

In principle, the fluctuations of Wr depend on the strength of the repulsive interactions, the 
number of particles N, the integration time r, the force ratio 9/a/ {STrVmax) and the depth Vmax 
of the substrate potential. Nevertheless, using the Langevin model of Eq. (H, we can gain some 
insight into the statistical properties of the work by analyzing two limit ideal cases which bear 
resemblance to static and sliding friction, respectively. 

The first case corresponds to N non-interacting particles moving under the influence a very 
weak force / <C SirYmax/i^o,) in presence of a very high potential barrier Vmax ^ ksT, such 
that the inverse Kramers rate of each particle becomes much larger than the other characteristic 
time-scales of the system. In this case, which resembles static friction conditions, the system is 
in a quasi-equilibrium state, where the particles are pinned by the potential wells at an average 
distance 9/a^/(167r^V)nax) from the minima in order to balance the external force /. The 


probability density function of the work w^ 
with mean fir and variance cr^ given by 


is Gaussian, i.e. P{wr) = 


\/2na4 


exp 


2^? 


Mr 


where k = [47r/(3a)]^V)nax is the effective stiffness of the restoring force exerted by a periodic 
light potential with hexagonal symmetry m- Note that, while fir = 0 because no mechanical 
work is done in average by f, is non-zero. This is due to the thermal fluctuations of the 
solvent molecules, which can promote either positive or negative work fluctuations by randomly 
moving the particles with or against the applied force. 

The second ideal case is when N non-interacting particles are driven by a sufficiently large 
force / 1:|> SnVmax /i^ci), such that they move at the highest possible average velocity (x) = f jj, 
where the influence of the periodic potential is negligible, similar to free sliding friction. In such 
a case, the probability density function of Wr is also Gaussian, with mean and variance 

fJjq- 


^ The linear relation fir oct is not necessarily fulfilled for systems with many coupled degrees of freedom driven 
by time dependent forces, see for example [37]. 


Nf 

ksTj 

2Nf 

ksT'j 


T = 2flr 


( 10 ) 
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Figure 4. Probability density function of the work done by a constant force / = 12.4 fN over 
T = 3.6 s on subsystems composed of different number n of particles across two potentials of 
depths (a) Vmax = SOksT and (b) Vmax = l^ksT. From top to bottom in|U[a) and from left to 
right in |l|^b); n = 1 (dark blue), 5 (light blue), 10 (cyan), 25 (dark green), 50 (light green), 100 
(orange), and 250 (red). Insets: work distribution rescaled according to Eq. (fT^ . The symbols 
correspond to n = 1 (+), 5 (*), 10 (o), 25 (x), 50 (□), 100 (v); and 250 (o). The solid lines 
represent the rescaled Gaussian distributions for non-interacting particles given by Eqs. ([9]) and 
m, respectively. 


respectively. We point out that only in this particular case, the non-equilibrium work trivially 
satisfies the detailed steady-state fluctuation theorem [2] 


P{wr = w) 

in —-- 

P{Wr = —W) 


2/7'r 

= - = W, 


( 11 ) 


because Wr is actually equal to the total entropy production of the system, normalized by ks- 
In the following, we discuss how the previous ideal expressions for P{wr) compare to the 
experimental work distributions for interacting particles in the static and sliding friction regimes. 
Furthermore, we also investigate the work distribution in the intermediate regime when tuning 
the control parameters to induce a transition from static to sliding friction of the colloidal 
monolayer. 
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Figure 5. (a) Probability density function of the work done by a constant force / = 12.4 fN on 
n = 50 particles against a periodic light potential of depth Vmax = over different time 

intervals: r = 0.6 s (+), 1.2 s (*), 2.4 s (•), 4.8 s (x), 9.6 s (□), 19.2 s (y), and 38.4 s (<). 
The solid line represents the case of non-interacting particles under the same 14iax and /. Inset: 
Standard deviation of Wr as a function of r for non-interacting (solid line) and magnetically 
coupled (dashed line) particles (b) Probability density function of the work done by a constant 
force / = 12.4 fN on n = 50 particles against a light potential of depth Vmax = over 

different time intervals. From left to right: r = 0.6 s, 1.2 s, 2.4 s, 4.8 s, 9.6 s, 19.2 s, and 38.4 s. 
Inset: work distribution measured over different r and rescaled according to Eq. (USD. Same 
symbols as in Fig. [5Ka). The solid line represents the case of non-interacting particles. 


5. Work distribution for interacting particles under applied force 

5.1. Subsystem size 

Since only a portion of the complete monolayer of ~ 5000 particles is accessible for data 
analysis, a frequent problem encountered in spatially extended systems [Ml Ea SO], we first 
investigate the effect of measuring the work done on a smaller subsystem composed oi n < N 
particles, thus ignoring its coupling with the N — n degrees of freedom of the rest of the system. 
A possible way to probe the role of such a coupling is by means of the differences between the 
statistical properties of the work applied on a subsystem of n interacting particles with those 
observed in a subsystem of the same size n of non-interacting ones, where there is no coupling. 
Note that in absence of interactions, the work is a Gaussian variable with mean and variance 
proportional to the number of components n of the subsystem for the two limit cases described by 

















Eqs. ([9]) and (fTOjl . i.e. the width of the distribution scales in both cases as ar oc ^/n. Therefore, 
upon translating the work distribution to the origin by an amount and then squeezing it by 
its width 


w 


* 

r 



P*{w*) = V^P{^W*+ ^Xr) 


( 12 ) 


any subsystem composed of n non-interacting particles exhibits a n-independent prohle P*{w*). 
This means that the statistical properties of the work done on the whole system can be probed 
by measuring the work done on any subsystem of arbitrary size. This situation can change 
drastically in presence of particle interactions, though. As discussed in US m], because of 
the spatio-temporal correlations created by the interactions between the subsystem and the 
surroundings, the undercount of slow degrees of freedom can give rise to strong modifications 
of the statistical properties of the thermodynamic quantities of the subsystem with respect to 
those of the complete system. Then, it is not expected that the variance scales as oc n for 
sufficiently small n in presence of interactions. The effect of the coupling with the surroundings 
only vanishes when the size of the sampling subsystem spans a length-scale larger than the typical 
correlation length induced by the interactions, thus recovering the actual statistical properties 
of the complete system [40] . Indeed, in presence of repulsive interactions we observe this kind 
of non-trivial dependence of the work distribution on the number n of NESS trajectories used 
to compute Wr from Eq. ([7|) in both static and free sliding frictional regimes. 

In EigHKa) we plot the probability density function P{wr) of the work computed over r = 3.6 s 
for subsystems composed of different number n of particles (n = 1,5,10,25,50,100,250) at 
/ = 12.4 fN and Vmax = SOksT, for which all the particles in the monolayer are pinned by the 
potential wells over the observation times accessible in the experiment, see left panel of Eig. [2](a). 
Each subsystem is chosen in such a way that n neighboring particles cover an approximately 
square area ~ na^. We find that for all the values of n, P{wr) is symmetric and peaked around 
w-r = 0 because no work is done in average in this quasi-equilibrium state, whereas its width 
increases with increasing n. In the inset of Eig. SKa) we plot the work distribution rescaled 
according to Eq. m- We observe that P*{w*) has a Gaussian profile, and unlike the case of 
non-interacting pinned particles, its width increases with increasing n. This implies that the 
variance ci^ grows faster than n in presence of repulsive interactions for sufficiently small n, an 
indication that the correlations between the subsystem and the N — n surrounding particles 
are signihcant. Nevertheless, for sufficiently large values of n, we hnd that P*{w*) seems to 
converge to a size-independent profile, as shown in the inset of Eig. S^a) for n > 100. The 
convergence demonstrates that finite-size effects due to the spatial correlations between the 
sampling subsystem and rest of the particles in the monolayer become negligible compared to 
the global behavior of Wr for sufficiently large n. However, the effect of the particle interactions 
on the fluctuations of w-r persists even for sufficiently large n. As a matter of fact, when 
comparing the experimental P*{w*) for n = 100 and 250 with that computed from Eq. Q with 
k = 2.2 X 10“®Nm“^ for non-interacting particles (solid line in the inset of Eig. UKa)), we find 
that the former are much wider than the latter. This suggest that, with increasing numbers 
of particles n, the randomness created by the strongly non-linear coupling accumulate, giving 
rise to fluctuations of w-r larger than those that would be otherwise observed in absence of 
interactions. 

A different behavior is observed for free sliding, where all the particles are able to move across 
the potential landscape at a mean velocity close to u = //y, as those shown in the right panel 
of Eig. [2](b). An example of such a behavior is shown in Eig. HKb) where we plot the probability 
density function P{wr) of the work done by a force / = 12.4 fN on subsystems formed by different 
number of particles, n = 1, 5,10, 25, 50,100, 250, across a potential of depth Vmax = 15A;bT. In 
this case, the work distribution is Gaussian, whose maximum is located at positive values of Wr, 
because the applied force is able to perform mechanical work by moving the monolayer. The 



mean work, which coincides with the location of the maximum of P{wr), increases linearly with 
increasing n, in quantitative agreement with Eq. ([8|), as shown in Fig. [3]^b). On the other hand, 
the presence of interactions affects the behavior of the fluctuations of Wr compared to the ideal 
sliding case described by Eq. (jlOp . In order to highlight these differences, in the inset of Fig. SH^b) 
we plot the rescaled work distribution P*{w*) defined in Eq. m- Once more, the effect of the 
correlation between the subsystem and the rest of the monolayer can be observed for small values 
of n, for which the width of the rescaled distribution increases with n. However, for n > 25, 
P*{w*) converges to a n-independent profile, thus probing the actual statistical properties of 
w-r for the complete system. This convergence implies that the variance of the work scales as 
a'^ (X n for sufficiently large n. In the inset of Fig. IHb) we also plot as a solid line the rescaled 
work distribution of non-interacting sliding particles, described by Eq. m- Interestingly, we 
find that the width of the rescaled work distribution in the presence of interactions is smaller 
than that of the non-interacting case. We can interpret this narrowing as a reduction of the 
work fluctuations due to the repulsive interactions, which give rise to an effective stiffening of 
the monolayer, thus preventing large random excursion of the particles induced by the thermal 
fluctuations around the drift imposed by f. This is consistent with the fact that for a perfectly 
stiff colloidal crystal, which can be realized in the limit of infinitely strong repulsive interactions, 
thermal fluctuations are suppressed [28], which gives rise to a complete sharpening of the work 
distribution around the mean value of Eq. ([8]). 

5.2. Integration time 

We now focus on the dependence of the probability density function of the work on the 
integration time r. We point out that for values of r smaller than the relaxation time-scales 
of the system, time-correlations can affect also the statistical properties of the work, because 
the expression of Wr in Eq. ([7|) involves differences at distinct times of the particle positions. 
Nonetheless, for sufficiently large values of r, such time-correlations vanish and therefore the 
shape of the work distribution must converge to a single profile upon time rescaling. 

We first show in Fig. [Sj^a) the results for the case of a pinned colloidal monolayer, where 
in average no mechanical work is done. Here we plot the probability density function of the 
work Wr done on n = 50 particles by a force / = 12.4 fN against a light potential of depth 
Vmax = SOksT over different integration times, 0.6 s < r < 38.4 s. We observe that, for all 
the values of r. P{wr) is Gaussian and centered around Wr = 0, whose width increases with 
increasing r. However, for r > 9.6 s, the width of the distribution levels off and all curves 
collapse onto a master curve regarless of r. This is further verified in the inset of Fig. [5]^a), 
where we plot as a dashed line the dependence of the standard deviation Ur of the work on r, 
observing a saturation to a constant value at sufficiently large r. The dependence of ar on r 
is qualitatively similar to that for non-interacting particles, shown as a solid line in the inset of 
Fig. [5]^a). In this case, according to Eq. Q, the variance of the Gaussian work distribution levels 
off exponentially for integration times larger than the viscous relaxation time of the particles in 
the potential wells, 'jjk = 2.8 s. This behavior of Wr can be actually understood at the single¬ 
particle level. For r <C 'j/k, the particle motion is strongly auto-correlated in time due to the 
energy stored by the confining light potential, which translates into a very narrow distribution 
P(wr). The motion becomes less and less correlated when r approaches "j/k, and therefore 
each particle is able to perform larger Brownian displacements within the potential wells both 
with and against the applied force, thus resulting in a broadening of P{wr). Nevertheless, the 
fluctuations of Wr cannot grow indefinitely with increasing r because the single-particle motion 
is always bounded to the potential wells, giving rise to a saturation of Ur for r > y/A: . Although 
qualitatively similar as a function of r, we observe a quantitative difference at r y/A: between 
the standard deviation of Wr for interacting particles with respect to that in the non-interacting 
case, as shown in the inset Fig. [5]^ a). This difference is due to the strong coupling between the 


particles, which gives rise to a complex non-linear particle dynamics within the potential wells. 

In Fig. [5]^b) we illustrate the dependence of the work distribution P{wr) on the integration 
time r for the sliding friction regime of n = 50 particles driven at / = 12.4 fN and Vmax = 

We find that P{wr) is Gaussian and the location of the maximum increases linearly with 
increasing r in accordance with Eq. ([8|). Once more, inspired by the comparison with non¬ 
interacting particles, where the width of the distribution scales as ar oc y/r (see Eq. (fTn|) l. we 
can test a scaling with respect to r similar to Eq. m 
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P*{w*) = y/rPiy/rW* + Hr)- 


(13) 


Interestingly, in the inset of Fig. EKb) we show that the work distributions, rescaled according 
to Eq. (I13p . collapse onto a master curve for all r. This essentially means that in this frictional 
regime the variance of the work scales as cr^ oc r even in presence of particle interactions. Note 
that in this case, the particles are not confined to move in the potential wells, and consequently 
there is no intrinsic relaxation time in the dynamics, which explains the very fast convergence 
of P*{w*) to the master curve. Quantitative differences are observed between the experimental 
P*{w*) and the case without interactions (solid line in the inset of Fig. [5Kb)), though. This occurs 
due to the narrowing of the work distribution due to the effective stiffening of the monolayer. 


5.3. Depinning transition 

We now show how the shape of the work distribution changes between the two very distinct 
cases previously studied, i.e. when changing the experimental parameters to induce a transition 
from the regime where all the particles are pinned on the substrate, to the depinning of the 
colloidal monolayer and subsequent free sliding. We recall that in the two extreme regimes of 
static and sliding friction, the work distribution is Gaussian, even though the mean and variance 
behave differently as a function of n and r. While in static friction these quantities also depend 
strongly on both the elastic stiffness k exerted by the substrate potential and the repulsive 
pair-interactions, they are only affected by the strength of the interactions for sliding friction. 

In Fig. EKa) we illustrate the effect on the shape of the work distribution for n = 50 particles, 
computed over r = 9.6 s, when increasing the value of the applied force / at constant potential 
depth Vmax = SOksT in order to induce particle depinning. Interestingly, we observe that the 
work distribution becomes asymmetric with respect to the maximum with increasing /, as can 
be observed for / = 18.6 fN and / = 24.8 fN. In particular, non-Gaussian tails appear at positive 
values of Wr, as highlighted in the semilog plot of the inset of Fig. (IKa). For these values of /, 
the spatio-temporal dynamics of the monolayer becomes heterogeneous, as can be observed from 
the particle trajectories of Fig. E^a). For example, for / = 18.6 fN there are regions where the 
particles are still confined by the potential wells, because the external force is still smaller than 
the maximum conservative force exerted by the light field: 9fa/{STrVmax) = 0.54. However, 
the combination of thermal fluctuations and non-linear repulsive interactions can promote hops 
to the neighboring potential wells, thus creating collective motion of clusters of particles. The 
collective motion is in turn facilitated by the symmetry beaking induced by the external force. 
While the stagnant particles do not contribute to the mean value of the work but only to the 
fluctuations around Wr = 0, the non-Gaussian tails originate from the work done on the sliding 
particles. This regime persists even when the complete monolayer can slide across the periodic 
substrate, as observed at / = 24.8 fN and Vmax = SOksT, where the particles undergo stick-slip 
motion. In this case, the spatial heterogeneity is induced by the corrugation potential, creating 
zones around the potential minima where the particles slow down, whereas they move faster 
when overcoming the potential barriers, as illustrated by the trajectories in the right panel of 
Fig. EKa). Note that this heterogeneity results also in a mean particle velocity (x) = 105nms“^ 



much smaller than the maximum velocity that could be achieved in presence of a completely 
flat substrate {v = 400nms“^). 

Gaussianity of Wr is recovered at sufficiently large 9fa/{87rVmax), though, with a respective 
narrowing of P(wr)- This can be observed in Fig. EJ^b), where we plot P(wr) for the second 
experimental protocol with which we can reach more easily the free sliding regime. A reduction of 
only in the potential depth Vmax, from 29k bT to Ih/csT, is enough to observe a prominant 

change of the shape of the work distribution, as shown in the inset of Fig. El^b). For this values of 
Vmaxi the mean particle velocity changes from 74nms“^ to 185nms“^, whereas the maximum 
velocity that could be achieved for this value of / on a completely flat surface is 200 nm s“^. This 
dramatic change in the shape of P{wr) reveals that not only the average tribological properties 
of the monolayer [26| but also the properties of the fluctuations of the work done on it become 
very sensitive when tuning the experimental parameters close to the depinning transition. 


6. Asymmetry of the non-equilibrium work distribution 

Finally, we investigate the work distribution P{wr) within the context of non-equilibrium work 
relations. More specifically, we focus on the evaluation of quantity In , which quantihes 

the asymmetry of the probability of observing positive work fluctuations, where the monolayer 
moves in the direction of the applied force, with respect to the probability of observing rare 
negative fluctuations, where the monolayer moves against the force. We point out that, although 
empirically satisified in many steady-state complex systems [Ml sa Ha ig sa sg El Hsi sg EH] , 
a simple linear relation such as Eq. m for the asymmetry function. In oc w, is not 

expected to hold generally for the NESS system we study. Indeed, for this kind of frictional 
processes, such a linear asymmetry relation is strictly valid only for the work done by an external 
force on a collection of non-interacting particles freely sliding on a perfectly flat substrate. In 
such a case, Eq. m is a direct consequence of the detailed Fluctuation Theorem, which only 
applies to the total entropy production of the system, and which in that specific case equals 
the work done by an external force. On the other hand, since the work defined by Eq. © 
has a definite parity under time-reversal, it must satisfy a generalized Fluctuation Theorem in 
presence of the corrugation potential, the pair-interactions and the global confining potential of 
the sample cell |2] 


In 


PjWr 

P{Wt 
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(14) 


In Eq. (fT4|l . Acr = [E{t) — E(0)]/(/cbT) is the variation of the total potential energy of the 
system during a time interval r, given by Eq. m, Asr is the stochastic entropy change over 
r m, and the brackets denote a conditional average over the stochastic realizations for which 
Wr equals the value w. The last term on the right-hand side of Eq. m is in general non-zero 
in presence of particle interactions, and therefore P{wr) does not necessarily satisfy the exact 
linear relation of Eq. m- Note that if P{wr) is Gaussian, the asymmetry function can still be 
proportional to w, i.e. 


, P(Wr = +W) 

In —--- 

P{Wr = —Wj 


aw, 


(15) 


where the parameter a, i.e. the slope of the linear relation, is given by 



(16) 


However, unlike the ideal case of non-interacting particle described by Eq. m, the parameter 
a is in general different from 1 because the second term on the right-hand side of Eq. 
which involves the particle interactions and the substrate potential, is non-zero. Eor instance. 








[units of k^T] 


Figure 6. (a) Probability density function of the work done by different external forces / on 
n = 50 particles against a potential of depth Vmax = SO/c^T, computed over r = 9.6 s. From 
left to right: / = 12.4 fN, 18.6 fN, and 24.8 fN. (b) Probability density function of the work 
done by an external force / = 12.4 fN on n = 50 particles against potentials of different depth 
Vmax, computed over r = 9.6 s. From left to right: Vmax = S0kBT,25kBT,20kBT, and Ih/csT. 
The insets are semi-logarithmic representations of the same plots in the main figures. 


we observe that Eq. (I15p holds in the static friction regime (completely pinned monolayer) and in 
the sliding regime, as shown in Fig. [7Ka). In the static friction regime, the parameter a is equal 
to 0 for all n and r, because the system is in a quasi-equilibrium state, with equal probabilities 
P{wr = +w) and P{wr = —w). On the other hand, we find a ~ 2 for all r in the free sliding 
regime, as illustrated by the symbols around the dashed line in Fig. [TKa). This implies that in 
this case the second term on the right hand side of Eq. (|14p is non-zero: ~ —w. It 

should be noted that the direct computation of the asymmetry function from P{wr) is restricted 
to rather small values of n and r, because negative work fluctuations are difficult to sample with 
increasing values of such parameters. However, taking into account that P{wr) is Gaussian, ol 
can be estimated from the mean and the variance by means of Eq. m- Surprisingly, in the 
inset of Fig. EKa) we show that the value a ~ 2 holds even for n and r as large as 500 and 40 s, 
respectively, thus demonstrating that the particle interactions give rise to a robust behavior of 






















the term ~ —w in Eq. ()14h . This unconventional behavior of the asymmetry 

function can be traced back to the strong coupling between the particles forming the crystalline 
monolayer. Indeed, with increasing strength of repulsive interactions, which is fixed in our 
experiment by the magnetic field B, the width of the work distribution decreases because of the 
increasing stiffness of the colloidal crystal. Note that in the limit of a perfectly stiff colloidal 
crystal, i.e. created by infinitely large repulsive interactions, the fluctuations of w-r are completely 
suppressed. In this case the work distribution becomes a delta function, P{wr) = 6{wr — Hr), 
which gives rise to a —)■ oo. Hence, for sliding friction a must be an increasing function of the 
pair interaction strength, bounded by the values a = 1 (no interactions) and oo (infinitely large 
repulsions). The value a ^ 2, specific to our experimental conditions, clearly illustrates that 
repulsive pair interactions reduce the fluctuations of Wr compared to the value a = 1 of Eq. m 
in absence of interactions. 

The intermediate frictional regime, where P{wr) exhibits non-Gaussian tails due to the 
heterogeneous spatio-temporal dynamics of the monolayer, is particularly interesting. In this 
case, we observe that the asymmetry function is not even linear in w, as shown in Eig. [T^b) for 
/ = 12.4 fN and Emax = 25kBT, at which only partial depinning from the substrate is achieved. 
The asymmetry function is approximately linear for small values of w, with slope a <C 1 due to 
large negative work fluctuations on the stagnant particles. Nevertheless, significant deviations 
from this linear behavior show up at larger work fluctuations, w > 5, when probing values of 
Wr on the non-Gaussian tails plotted in Eig. [IKb). Once more, the behavior of the asymmetry 
functions seems to be robust, as shown in Eig. EKb), where all the data points collapse to a master 
curve for different values of r and n. We point out that in general, such a non-linear behavior of 
the asymmetry function is not easily observed in systems described by a small number degrees 
of freedom, because in such a case large negative fluctuations are difficult to sample [H]. In our 
experimental system we are able to achieve this because of the existence of strong negative work 
fluctuations, which originate from the heterogenous dynamics of the coupled degrees of freedom 
of the system during the depinning transition of the colloidal monolayer. 

7. Summary and conclusion 

We have investigated the statistical properties of the work done by an external force on a 
monolayer of magnetically interacting particles driven across a periodic potential, which mimics 
friction between crystalline atomic surfaces. We have studied the influence of the number of 
particles used to probe these properties, the integration time, and the control parameters that 
are tuned to induce a transition from a pinned state (static friction) to complete depinning 
from the substrate potential (sliding friction). We have shown that, in the static and free 
sliding regimes, the work distribution converges to a Gaussian master curve for sufficiently 
large number of particles and integration times upon rescaling of these parameters. We have 
found that the mean and variance of such work distributions depend on the strength of the 
repulsive interactions, which in particular give rise to a stiffening of the monolayer for free sliding. 
Interestingly, we have also found that in the intermediate friction regime, where the monolayer 
undergoes a depinning transition, the work distribution becomes non-Gaussian because of the 
heterogeneity of the particle dynamics, e.g. due to partial depinning and stick-slip motion. We 
have shown that in general, the work distribution exhibits unconventional asymmetry properties 
within the context of non-equilibrium fluctuations relations. We have demonstrated that such 
a behavior originates from the presence of repulsive particle interactions. Thus, we provide 
the first experimental measurements of a stochastic thermodynamic quantity with non-trivial 
properties for a mesoscopic system with many coupled degrees of freedom. 
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Figure 7. (a) Asymmetry function of the work distribution P{wr) for: static (+), and sliding 
friction regime measured for n = 10 particles at / = 12.4 fN, Vmax = over r = 0.6 s 

(o), r = 1.2 s (□), T = 2.4 s (o). The solid circles are measurements under the same conditions 
for n = 25 particles during r = 1.2 s. The dashed line is a guide to the eye with slope a = 2. 
Inset: dependence of the parameter a on r computed by means of Eq. (1161) for n = 100 (i>) and 
500 (solid line) particles, (b) Asymmetry function of the work distribution P{wr) measured for 
n = 25 particles moving with a heterogeneous dynamics under / = 12.4 fN and I4iax = 25 A:bT 
over r = 0.6 s (+), r = 1.2 s (*), r = 2.4 s (o), r = 4.8 s (x), r = 9.6 s (□), r = 19.2 s (o), 
and r = 38.4 s (o). The solid circles are measurements under the same conditions for re = 50 
particles over r = 38.4 s. The dashed line is a guide to the eye with slope a = 0.15. 
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